Geometrical frustration induced (semi-) metal to insulator transition 
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We study the low-energy properties of the geometricaUy frustrated Ifubbard model on a three- 
dimensional pyrochlore lattice and a two-dimensional checkerboard lattice on the basis of the renor- 
malization group method and mean field analysis. It is found that in the half-filling case, a (semi- 
)metal to insulator transition (MIT) occurs. Also, in the insulating phase, which has a spin gap, 
the spin rotational symmetry is not broken, while charge ordering exists. The results are applied to 
the description of the MIT observed in the pyrochlore system TI2RU2O7. 



Recently, the role of geometrical frustration in both 
localized and itinerant electron systems has attracted re- 
newed interest |0,||j3|J|,^,^,0j8) . In localized systems, the 
possibility of exotic phase such as a spin liquid with- 
out long-range magnetic order has been extensively in- 
vestigated both experimentally ||l|,||j^J^ and theoretically 
l^'iHlJll- itinerant systems, the presence of charge de- 
grees of freedom provides a route for the relaxation of 
magnetic frustration. However, when electron correla- 
tion is sufficiently strong, the magnetic frustration may 
still affect the low-energy properties significantly. For 
example, it has been found experimentally that some py- 
rochlore oxides, such as TI2RU2O7 and Cd20s207, ex- 
hibit a (semi-)metal to insulator transition (MIT) at fi- 
nite critical temperatures |^,^. Since such systems pos- 
sess the fully-frustrated lattice structure, referred to as 
a network of corner-sharing tetrahedra (that is, a py- 
rochlore lattice), the magnetic properties of the insu- 
lating phase are yet largely unexplained. Moreover the 
mechanism of the MITs observed in these systems is still 
an open problem. Geometrical frustration may play an 
important role in the MITs. From this point of view, in 
the present paper, we study the interplay between elec- 
tron correlation and geometrical frustration in the Hub- 
bard model on a three-dimensional (3D) pyrochlore lat- 
tice and on a two-dimensional (2D) checkerboard lattice, 
the so-called 2D pyrochlore (FIG.l). Although real py- 
rochlore oxides have electronic structure composed of t2g 
orbitals, the present study on these simpler single-band 
models may provide important insight into the role of 
geometrical frustration in MIT. Furthermore, TI2RU2O7 
has, apart from the t2g band, a nearly half-filled Tl s 
band, whose important features are described by the 
3D pyrochlore Hubbard model We believe that this 
model may provide a useful understanding of the MIT 
undergone in this material. 

The non-interacting energy bands of these two Hub- 
bard models have a common interesting feature: They 
consist of a flat band (or two degenerate flat bands) 
on the upper band edge and a dispersive band that is 
tangent to the flat band (or flat bands) at the F point 



|iO| . According to graph theory, this band structure 
stems from the geometrical frustration inherent in the 
pyrochlore lattice. In the half-filling case, n = 1, on 
which we concentrate henceforth, the dispersive band is 
fully occupied, and the flat band(s) is empty. In the non- 
interacting half-filling case, the system is in a semi-metal 
state, since the Fermi velocity is vanishing, though there 
is no excitation gap. We study how this state is affected 
by electron correlation. 
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FIG. 1. 2D and 3D pyrochlore lattices. 

Diagonalizing the kinetic term, we write the Hamilto- 
nian as 
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where m = 2 in the 2D case, and to = 4 in the 3D 
case. In the 2D case, i?fci = 2, E-k2 — 4 cos fc^; cos fcj, — 2, 



s\\(k) — S22(fc) = sin( "'^ " )/-^l — cosfc^ cosfc^, and 
S2i(fc) — — si2(fc) — sin( '''' ) / -^1 — cos cos ky . In 



the 3D case, E^^ = Ek2 - 2, Eks = -2 + 2^/1+1^, 
and Ek4 — —2 — 2yT+^, with tk — cos(2/cj;) cos{2ky) + 
cos(2fcj,) cos(2fc2)-t-cos(2fcz) cos(2fca;). The form of s^^{k) 
in the 3D case is given in Ref 0. The annihilation op- 
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erator of electrons at the /i-th site in a unit cell is given 

by Ckfj,a = J27=i Sij.v{k)akva- 

As shown in Ref the perturbative calculation in 
U for the above Hamiltonian suffers from divergences at 
third and higher order, due to the presence of the flat 
band(s). To treat the divergences in a controlled man- 
ner, we apply the renormalization group (RG) method. 
In previous studies of electron systems pl| , a momentum 
cutoff that separates the neighborhood of the Fermi sur- 
face from the higher momentum part is introduced. How- 
ever, in the presence of the flat band(s), this procedure is 
not applicable. To overcome this problem, we introduce 
the infrared energy cutoff A in the following manner: 
''Pf.aik.en) = ^>Ak,en)Q{\en\ - A) -h V^^(fc,e„)e(A - 
|e„|). Here, 'ip^a{k,e) is the Grassmann field correspond- 
ing to flfepcr- 

Using a standard method, we can obtain the RG equa- 
tions of the single-particle self-energy, S^i,(fc), and the 
4-point vertex functions, rQ,^^5(fci, ^2; ^3, ^4) up to the 
one- loop order. In our systems, there are six species 
of 4-point vertices, as shown in FIG. 2(a), apart from 
the spin degrees of freedom and the two-fold degen- 
eracy of the flat bands in the 3D case. We assume 
that the momentum dependences of the 4-point vertex 
functions are given mainly by r°(fci, ^2; fca, /C4) in the 
renormalization processes. This is made explicitly by 
replacing ra&afc(fci, ^2; ^3, ^4) with gir[Jf,^j(fci, fe; ^3, ^4), 
Ybbba{ki,k2\k'i,kA) with g4rgf,ba(^i' ^2; ^3, ^4), etc. This 
approximation is fairly good, because in the vicinity of 
the r point, where the most important scattering pro- 
cesses occur, the band structure is almost isotropic. Be- 
cause the flat bands are empty and the dispersive band is 
fully occupied, the particle-particle processes between the 
flat bands and the particle-hole processes between the flat 
bands and the dispersive band give the leading singular 
contributions (see FIG. 2(b)). We take into account these 
contributions in the derivation of the RG equations. We 
found from the analysis of the RG equations that, among 
the six running couplings, 52, 53, and 55 are irrelevant in 
the low energy limit. The RG equations for the other 
couplings, f/i, 34, and f/g are written 
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where r/ = (d - 2)/2, / = ln(Ao/A), with Aq the band 
width, and d is the spatial dimension. The couplings gis 
and git denote the spin singlet and triplet parts, respec- 
tively. In the 2D case, a = Yl,ki^ii{k) — Sii(fc)si2(^))/2, 
h = l/(32t), and in the 3D case, a = J2kli4iik) + 



s?2(fc))' - (sii(fc)s2i(fc) + si2(fc)s22(fc))2]/2, and b « 
In the derivation of these equations for the 
3D case, we have used the fact that in the vicinity of the 
r point, the two degenerate flat bands do not mix with 
each other in scattering processes. Thus in this case the 
two-fold degeneracy just gives an overall factor of 2. 
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FIG. 2(a) The six species of 4-point vertices. Here "a" 
and "6" indicate the dispersive band and the flat band, 
respectively, (b) The leading singular bubble diagrams. 

2D checkerboard pyrochlore- We first consider the 2D 
case, whose theoretical treatment is simpler. We solved 
the RG equations (^-(^ numerically for a particular set 
of parameter values, and obtained the RG flow shown 
in FIG. 3(a). We found that for any small value of U/t, 
git flows into the strong-coupling regime. This indicates 
some instability in this channel. Although gis also scales 
into the strong-coupling regime, it is sub-dominant com- 
pared to git- We also show in FIG. 3(a) the RG flows of 
the couplings 3git — gis and gis + gu, which are related 
to the charge and spin susceptibilities, respectively. We 
see that some instability appears in the charge degrees of 
freedom. To elucidate the nature of this instability, we 
write the RG equations for the single-particle self-energy 
in the following form: 
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26(3.gu-5i.)(S^2TT+Smi) 
26(gis + git)(E^2TT ~ ^]^2i i) 
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dl 



2&(.9is +5it)S]^2n- 



(6) 
(7) 
(8) 



In the derivation of these equations, we have ignored the 
diagonal self-energy, which are not important in the fol- 
lowing argument, and expanded the RG equations up to 
the flrst order in S^2- Because the strongest divergence 
of the 4-point vertex appears in 3git — gis (see FIG. 3(a)), 
the off-diagonal self-energy ^i2crcr becomes non-zero 
at some critical A = Ac. This is easily seen by solving (H), 
which gives E.S^2.. - E.sXexp[26/o'dr(35it - 
gis)]- Note that ^i2a-a vanishing in the vicinity of the 
r point, because of the momentum dependence of Sp^uik)- 
Thus for A = Ac at which value 3git — gis is divergent, 
J2a '^i2a-a- bccomcs nou-zcro. 

The above RG analysis implies the existence of a mean 
field solution for which the order parameter is given by 



2 



yO,k2cj)- This state is 



characterized by electron-hole pairing with parallel spins, 
which leads to the formation of both spin and charge 
gaps preserving the spin rotational symmetry. Accord- 
ing to the numerical analysis of the RG equations (||)-(||), 
gus is mainly renormalized by the first term of the right- 
hand side of (|^). Then, the renormalized coupling g^s is 
approximately given by RPA-like expressions. As a re- 
sult, we obtain the self-consistent gap equation for Afc 
expressed diagramatically in FIG. 4. The transition tem- 
perature is determined from the linearized gap equation, 

Afc = ;^n(fc,(7-fc)G„((7-fc)G22(<z-fc) 

q,k' 

X n(fc',g~fc)Gii(fc')G22(fc')Afc', (9) 

where G^p(fc) = l/(e„ - Ek^j), k = (i£,fc), and 
n(fc,fc') = Y.u=±vUt''{k,k')/{l - c,D{k + fc')), 
t±(fc,fe') = (sn(fc)si2(fc') ± si2(fc)sn(fc'))V2 with 
D{q) - -T;7E„,fcGii(fc)Gii(g-fc), c+ = 2a, and c_ = 
'^kSii- Here we have ignored the diagonal self-energy. 
Equation implies that the gap function can be written 
Afc = sii(fc)si2(fc)Ao, where Aq is a constant. From (H), 
we have Aq = AQlu^/m)[\n{8t/U)-bo] lii{8t/Tc), where 
bo = 0.322. For U < Uc 0.725 x 8t, a state with non- 
zero Aq is realized. We have also applied the Ginzburg- 
Landau analysis to this mean field solution and found 
that, in the 2D case, the transition temperature vanishes 
in accordance with the Mermin-Wagner-Coleman theo- 
rem. Nevertheless, the above analysis demonstrates that 
in the ground state at zero temperature, the gap Afc is 
non-zero, and the system is in an insulating state. 

We now further investigate the properties of the in- 
sulating phase. In this phase, because the order pa- 
rameter does not break the spin rotational symmetry, 
there is no long-range magnetic order. However, a spin 
gap exists. The spin-spin correlation function obtained 
from the above mean field solution is Imxs(9,w)/w ~ 
(l/(2rcosh2(Afc/2r)))fc.^o for uj ^ 0. Here, (• • •)fc^o is 
the angular average near fc = 0. This spin gap behavior 
can be observed using NMR measurements. 

Another important property of the insulating phase in- 
volves the charge degrees of freedom. The formation of 
the gap Afc gives rise to a difference between the charge 
densities at the sites 1 and 2 in a unit cell given by 
pi — P2 ^ Ao/t, up to a constant factor. Thus charge or- 
dering (CO) with a charge density displacement propor- 
tional to the gap characterizes this insulating state (see 
FIG. 5 (a)). This noteworthy result can be understood as 
follows. In our system, three electrons occupying nearest 
neighbor sites cost energy loss caused by magnetic frus- 
tration. Conversely, magnetic frustration induces an ef- 
fective finite-range repulsion between electrons at nearest 
neighbor sites. If this finite-range repulsion is sufficiently 
strong to overcomes the on-site Coulomb interaction U, 
the CO state will be stabilized. This is possible if U 



is not so large. As U increases, a transition to a con- 
ventional Mott insulator with no charge ordering should 
occur. This transition cannot be described within our 
weak-coupling analysis. 
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FIG. 3(a) The RG flow of the running couplings in the 
2D case with U/8t = 0.25. (b) The RG flow in the 3D 
case with U/8t = 0.75. 
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FIG. 4 Diagrams for the linearized gap equations. 

3D pyrochlore- The above analysis can be straightfor- 
wardly generalized to the case of a 3D pyrochlore lat- 
tice. We obtain the RG flow numerically from (|^)-(|^) for 
d = 3. Here, in contrast to the 2D case, for sufficiently 
small U all couplings are irrelevant, and thus the semi- 
metal state is stable. However, for values of U larger 
than a certain critical value but still smaller than the 
band width, RG flow similar to that in the 2D case is 
obtained, as shown in FIG. 3(b). The coupling Sgu — 5is, 
which is related to the charge degrees of freedom, scales 
into the strong-coupling regime. This RG flow implies 
that, as in the 2D case, a particle-hole pairing state with 



order parameters A^ 



(13) 



^fclcr 



akscr) and A 



(23) 



'^ai^k2cr'^>^3'^) realized. Although the value of U used 
here is relatively large, we expect that the one-loop RG 
calculation still gives qualitatively correct results, as long 
as U is smaller than the band width. To examine the 
validity of the one-loop calculation, we explore the self- 
consistent mean field solution. The self-consistent gap 
equations for the 3D case are also given by the dia- 
gram shown in FIG. 4, from which we find that the gap 



functions are given by A 



(13) 



(13) 
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(23) 



Using the symme- 
try properties of s^v{k) in momentum space, we can 
show without solving the gap equations that Al, = 



0, and A\ 



(23) 



= A 



(23) 



= A 



(23) 



Thus, from the or- 



thogonal relations among the Sp,y(fc), we have A^^"^-* — 



(23) N 



The quantity Af 



Ai- 



(23) 



S42(fc)s43(fc)(Af^-A 

is determined from the gap equation. According to the 
RG analysis, the transition occurs only for sufficiently 
large U. Therefore, to determine the transition tempera- 
ture and the gap function correctly, we need to take into 
account the self-energy corrections i.e., pair breaking ef- 
fect. This calculation is rather involved, and we have not 
yet carried it out. However, we see from the RG flow that 
at the critical temperature Tc ^ Aoe~''= — 8t x 0.004, a 
transition from a semi-metal to an insulator occurs. In 
the resulting insulating state, the three-fold degeneracy 
at the r point in the semi-metal state is lifted completely, 
and a spin gap as well as a charge gap exists. 



yet been investigated experimentally. The experimental 
determination of whether a GO state exists for TI2RU2O7 
is a crucial test of this theory. When there exists coupling 
to a lattice, CO should accompany lattice distortion. It 
has been found experimentally that, in TI2RU2O7, the 
lattice structure changes from cubic to orthorhombic at 
the MIT point. This observation seems to suggest the 
presence of large charge fluctuations in this system. 

In conclusion, the 2D and 3D pyrochlore Hubbard 
models at the half-filling show the transition from semi- 
metal to spin-gapped insulator. In the insulating state, 
charge ordering occurs so as to relax geometrical frustra- 
tion. 

The author is grateful to K. Yamada and H. Tsunet- 
sugu for comments and discussions. This work was sup- 
ported by a Grant-in- Aid from the Ministry of Education, 
Science, and Culture, Japan. 




FIG. 5 (a) The CO pattern in the 2D case, (b) The GO 
pattern in the 3D case. 

As in the 2D case, we examine now the possibility of 
a GO state in the 3D system. In this case, there are four 
sites in a unit cell. The appearance of a gap causes a 

charge density displacement on each site given by 6pi, = 

^ (23) 

2Z]fc 5i'2(fc)si/3(fc)A^ for ^ = 1,2,3. Using the symme- 
try properties of s^u{k), we find 6pi = Sp2 = Sp^ ^ and 
Sp4 = —ZSpi. It is thus found that CO with the pattern 
displayed in FIG. 5(b) occurs in the insulating phase. In- 
terestingly, a similar CO pattern is observed in the spinel 
system AIV2O4 which possesses a V-site corner-sharing 
tetrahedron network |p^ . 

We now apply the results obtained above to the de- 
scription of the MIT observed in TI2RU2O7. The band 
calculation gives the band width of this system 8i ~ 2eV 
1^. Experimental data on the size of U do not exists. 
However, typically, the value of U for transition metal 
oxides is ~ 2eV. This gives us reason to believe that 
the analysis given in this paper, which suggests that the 
MIT occurs for large U ^ 8t, can be applied to the de- 
scription of the TI2RU2O7 system. The transition tem- 
perature estimated from the RG analysis is Tc ^ 96K, 
which is almost comparable with the experimental val- 
ues 100 '-^ 120K [Q. A recent NMR measurement has re- 
vealed the presence of a spin gap in the insulating state, 
which is consistent with our results ||l^. The possible 
existence of a CO state and large enhancement of charge 
fluctuations above Tc predicted in our theory have not 
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